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We explore the effect of mutual gravitational interaction between ultra-cold gas atoms on the 
dynamics of Bose-Einstein condensates (BEC). Small amplitude oscillation of BEC is studied by 
applying variational technique to reduce the Gross-Pitaevskii equation, with gravity included, to 
the equation of motion of a particle moving in a potential. According to our analysis, if the s-wave 
scattering length can be tuned to zero using Feshbach resonance for future BEC with occupation 
numbers as high as « 10 20 , there exists a critical ground state occupation number above which the 
BEC is unstable, provided that its constituents interact with a 1/r 3 gravity at short scales. 


I. INTRODUCTION 

Gravity is the weakest of all forces. This is essentially 
due to the smallness of Newton’s gravitational constant 
(or, equivalently, largeness of Planck mass), measured 
on scales larger than tens of km*. However, to resolve 
issues pertaining to naturalness and hierarchy problems 
in the Standard Model of particle physics, it has been 
conjectured that if large extra dimensions exist, the ef¬ 
fective gravitational coupling strength can be larger at 
sub-mm scales^. With the advent of exciting precision 
experiments involving Bose-Einstein condensation of al¬ 
kali atoms and molecules at ultra-low temperatures^, it 
is but natural to study effects of enhanced gravity ensuing 
from large extra dimensions (LED) on such macroscopic 
quantum phenomena. 

In this context, Dimopoulos and Geraci have proposed 
an interesting experiment to probe gravity at sub-micron 
scale through measurements of relative phase evolution 
rates in Bose-Einstein condensates (BEC) prepared in 
coherent superposition of states localized at two distinct 
potential wells, both situated near a moving wall of al¬ 
ternating gold and silver metal objects that form a pe¬ 
riodic massive source of gravity^. Similarly, Sigurdsson 
has suggested measuring fringe shifts of interfering pair 
of BEC falling past a long and narrow cylindrical mass in 
order to estimate modified transverse gravitational accel¬ 
eration, provided that the LED sub-mm scale is in excess 
of 0.01 mm£. 

Interestingly enough, the typical separation between 
atoms in ultra-cold gases is only about few 100 nm. This 
induces one to explore effects of mutual gravitational 
interaction between individual atoms of a BEC on its 
quantum dynamics, and ask whether such weak but long 
range forces can lead to instabilities. In this paper, we 
carefully examine some aspects of these ideas using vari¬ 
ational method. 


II. GROSS-PITAEVSKII EQUATION AND 
LARGE EXTRA DIMENSIONS INDUCED 
GRAVITY 


For N identical bosons constituting a dilute BEC 
at temperature T ss 0°A', the many body wavefunc- 


tion ^(rl, 7'2,.., rjv) describing the condensate can be ex¬ 
pressed upto a good approximation (assuming that the 
bosons interact weakly with each other) as, 


N 

^(r!,r2,..,?w) = 

3 =1 


where tp(f) is the normalized ground state wavefunction 
for a single boson. As each boson, in this case, is approx¬ 
imately in the same state, i){r) acts as the condensate 
wavefunction. 

In the T = 0°A' mean field approximation, dynami¬ 
cal evolution of the condensate wavefunction if(r, t) (nor¬ 
malized to unity) is, to a good extent, governed by the 
Gross-Pitaevskii equation, 


III-— = 

dt 


~V 2 +V ext +N j V{r~u)\if(u,t)\ 2 d\ 


ip(r, t) 


(1) 

where m, V ext (f) and V(r) are the boson mass, the trap 
potential energy required to confine the BEC and the 
interaction potential energy between two bosons, respec¬ 
tively. 

For the present purpose, the interaction potential en¬ 
ergy V in eq.(l) is a combination of s-wave scattering po¬ 
tential and the inter-bosonic gravitational potential en¬ 
ergy V g , so that, 


4 T) ^ 

V{r — u ) = 711 S 3 (r — u) + V„{\r — u|) . (2) 

m 

where a is the s-wave scattering length. 

Substitution of eq.(2) in eq.(l) results in the standard 
Gross-Pitaevskii equation (GPE)£, 


dif 

= 


-^-V 2 +V ext + NgW(r,t)\ J 
2m 



u\)\tf>(u, t)\ 2 d 3 u 


i>(r,t ) 


( 3 ) 


where g = 4vh a . 

^ m 

It is interesting to note that the quantum dynamics 
of a BEC, comprised of ultra-cold bosonic atoms an¬ 
chored to a planar honeycomb optical lattice and inter¬ 
acting weakly to one another via a contact interaction 
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much like the first term of the RHS of eq.(2), is described 
by a nonlinear Dirac equation^. Furthermore, the pseu¬ 
dospin degrees of freedom associated in this case with 
the two inequivalent sites of the sublattice display half 
integral spin angular momentum features, stretching the 
graphene analogy farther, even though the system is a 
bosonic one£. 

The GPE of eq.(3) can easily be derived from the fol¬ 
lowing action by demanding it to be stationary under 
infinitesimal variations of ip and ip*, 

S = J dt j d 3 r C (4) 

where the Lagrangian density C is given by, 

£ = } + £ v **- v * + v '“* w 2+ 


+"y #| 4 + yl i>\ 2 J V g (\f— u\)\ip(u, t)\ 2 d 3 u (5) 

Now we come to the gravitational potential energy V g 
appearing in eq.(5). In the framework of LED grav¬ 
ity, the hierarchy problem of the Standard Model can 
be ameliorated if (a) there exists a fundamental energy 
scale M*c 2 (« 1-1000 TeV, orders of magnitude less than 
the Planck energy = 10 19 GeV) for all interactions, and 
(b) there are additional sub-mrn scale spatial dimensions, 
so that the perceived weakness of Newtonian gravity on 
large scales in (3+l)-dimensional space-time is due to 
the gravitational field lines spilling into the hidden spa¬ 
tial dimensions^. In this formalism, the gravitational 
potential energy V g (r) between two point masses nil and 
W 2 separated by a distance r is given by, 


III. VARIATIONAL METHOD, GRAVITY AND 
BEC OSCILLATION MODES 

Solving eq.(3) with V g given by eq.( 6 ) is a nontrivial 
task. Instead, we take recourse to a variational method 
developed to study stability and low energy excitations 
of BECi£~— . In this approach, the parameters of a trial 
wavefunction ip tr is obtained by demanding that the ac¬ 
tion is extremized by ipt r - Since attractive contact in¬ 
teractions (i.e. a < 0 ) is known to cause collapse of 
BE C 11 ’ 13 for sufficiently large N, stability analysis with 
gravitational interactions included is worth studying. 

For this purpose, we consider a spherically symmetric 
trap potential, 

Vext = T^rnw o r 2 ( 8 ) 

and choose a normalized trial wavefunctiontL, 

iptr(r,t) = A{t)exp(-r 2 /2cr 2 (t))exp(iB{t)r 2 ) (9) 

where A(t), a(t) and B(t) are amplitude, width and 
phase parameters, respectively, that need to be deter¬ 
mined from extremization of the action (eqs.(4) and (5)). 
As ipt r is normalized, A(t) and <j(t) are related by, 

\A(t)\ 2 = (7w(f ))- 3 (10) 

so that, 

A(t) = ( v / 7 T 0 '(t )) —3/2 exp (ry(f)) (11) 

where 7 (t) is a time dependent phase. Substitution of 
eqs.( 8 )-(ll) in eq.(5) and carrying out the spatial integral 
thereafter lead to the following Lagrangian, 

L — [ dPT £, = fl'y + Lt n t d~ — 7=—7 -h 
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where m p i = y 77 is the Planck mass corresponding 
to the standard Newton’s gravitational constant G and 
-R* (n) is the radius of the extra dimensional n-torus given 

by, 


R*(n) 


\ 2/n + 

mpi \ n 

M * / 27rM*c 


(7) 


for n =1,2,... According to eq.( 6 ), closer one probes 
stronger is the gravity on scales smaller than R*(n). In 
the next section, we examine its implications on low lying 
excitations of BEC. 


where the gravity term is, 

Lim = ^ J d 3 r[ip(f,t)\ 2 J V g (\f-u\)\ip(u,t)\ 2 d 3 u 

(13) 

Using eq.( 6 ) for V g the above integral can be evaluated 
analytically for n = 0 and n = 1 cases so that, 


where, 


CD 0 

for n = 0 

(14a) 
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for n = 1 

(146) 
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Oil = 


NR*hc( to \ 2 ^ 1 (2k + 1)!! 
2\[2 \m p i) ^fc!2 2fc (2fc + l) 


NR*hc / to \ 2 
2 \m p i) 


(15 6) 


Extremizating the action entails Euler-Lagrange equa¬ 
tions dL/dq.j) — ( dL/dqj) = 0, for j=l and 2, with 
qi = B, q 2 = cr and L given by eq.(12) ("f(t) is non- 
dynamical as it appears only as an additive total deriva¬ 
tive term in eq.(12)). The equations of motion are, 


B(t) 


m a 
2 h a 


and, 


+ . 2h 2 ~ h 2 gN 

TiB +- B 2 - --j- 7 ^— l - 

to 2 toct 4 \\/2'k 3 I 2 o 3 


1 2 

-mw 0 


where, 


(16a) 


-/n(o-) 

(166) 


fn(c r) = 


ao 

3tr 3 


for n = 0 


2 ai 

3a 1 


for 77 = 1 


(16c) 


By combining eqs.(16a) and (16b), one arrives at the rel¬ 
evant equation needed to study small amplitude oscilla¬ 
tions in an ultra-cold cloud of bosons, 


777<T = —TOUIqCT + 


gN 


mcr 3 2v / 27t 3 / 2 <t 4 


- 2af n (a) (17a) 


Employing the following dimensionless quantities^ that 
make use of the BEC ground state scale -Jfi/mw o, 


v = 


\Jh/mwo 


T = W()t, P = 


Na 


Ti" yjh/mwo 


gN h 2 - 

=> 0 = 3/ 2 = —^h mw 0 P , (176) 

2v27t 3 / 2 to to- 

along with eq.(16c) in eq.(17a), we obtain, 
d' 2 v 1 P 

- 1~2 = ~ v + ~3 + — + Fn(v) (18) 

for 77 = 0,1, where the dimensionless gravitational accel¬ 
erations have the forms, 

F ° (V) = -\/5 "fe) Ugn/mwo ) (Wa) 

and, F l( ») = -|w(^) (^) »- 3 (1») 


The RHS of eq.(18) corresponds to an effective potential 
$ n (u) (tz = 0,1) given by, 


$n(^) = 


1 


P 

H- ^+vFo(v) for 77 = 0 (19c) 

3v 3 





vF\(v) 
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for 


77 = 1 


(19d) 


In order to study small amplitude oscillation modes, one 
needs to find the minima of $, 1 ( 7 ;). So, from $' n (v) = 0, 
the task here boils down to determining the zeroes of the 
quintic polynomial, 

v 5 - v - P - v 4 F n (v) = 0 (20) 


To estimate numerically the real positive roots vq of 
eq.(20) and the excitation frequencies proportional to 
VKM, we make use of typical experimental length 
scales, 

\]h/mwo ~ 10 -4 cm; ( c/wo)(m/m p i) 2 w3x 10~ 26 cm; 


a « 10 6 cm; i?*(l) « 200 ^m; h/(mc) « 1.6x 10 16 cm , 

(2f) 

having in mind a BEC comprising of 133 Cs for which 
(■ m/m p i) 2 = 4.9 x 10 -34 . 

Since both P and —F n (v) increase with N with the 
latter being negligibly smaller by orders of magnitude due 
to the smallness of ( m/m p i) 2 inspite of the other factors 
(see eqs.(17b),(19a,b) and (21)), it is obvious that the 
s-wave scatterings completely swamp the gravitational 
corrections to the excitation frequencies. The oscillation 
modes of such a problem in the absence of gravity had 
already been studied by Perez-Garcia et alJ^ 

To circumvent the dominance of binary s-wave scatter¬ 
ing one may, along with augmenting N, invoke Feshbach 
resonance^—. This effect enables experimentalists to 
tune the scattering length a magnetically, and reduce it 
to zero. Hence, with a vanishing P, in the 77 = 0 case (i.e. 
pure Newtonian gravity), one finds that for N < 10 21 , 
the real positive root vq of eq.(20) is very close to unity 
corresponding to a frequency of w = 2wq, as though the 
presence of F 0 (v) did not matter. 

However, for macroscopically large occupation num¬ 
bers N = 10 22 and 10 23 (BECs of future), one finds sig¬ 
nificant departures: vq = 0.78, uj = 2.4wo and vo = 0.12, 
u) = 66u>o, respectively. Because of the 3/vq term in 
4 >"(7 ;o), one expects a higher excitation frequency as vq 
becomes smaller than unity. Although these results sug¬ 
gest that rise in self-gravity due to increase in N beyond 
10 22 makes the ultra-cold gas cloud shrink drastically, 
caution needs to be exercised in concluding so. For, when 
the number density ~ N(\jTi/mwo vo)~ 3 becomes very 
large, other subatomic effects will start dominating and, 
also, it is likely that the variational method demands 
more care in such circumstances. For instance, when 
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N = 10 22 , our result vo = 0.78 implies a mean separa¬ 
tion between atoms in the condensate to be about 10 -11 
cm! Nevertheless, the observed pathology for N > 10 22 
situation suggests that it would be interesting to study 
the numerical solutions of GPE, with Newtonian gravity 
added, for macroscopic BEC. 

There is another way of getting around the problem 
of high density for large occupation numbers. One could 
increase the length scale y/Ti / (mwo) by choosing a weaker 
trap potential. Hence, to ensure mean separation not to 
fall below 10 Angstroms, the trap frequency Wq must 
satisfy the condition, 


Wq < 10 14 Uq 



h/m \ 
cgs units ) 


where vq is the positive root of eq.(20) corresponding to 
the Newtonian gravity case. 

In the n=l case (1/r 3 gravity), when P = 0, the non¬ 
zero roots of eq.(20) satisfy, 

= 1 — ( 2N/3)(m/m v i) 2 (R*/(h/mc )) (22) 

implying that the roots are complex when, 

N > N cr = (3/2){m/m p i)- 2 {R*/(h/mc))~ 1 (23) 

This is easily understood given that the potential of 
eq.(19d) can be expressed as, 


_ , , 1 2 1 
*■(■)= r + s? 


1 - 


N 

Wr 


(24) 


This signals instability for the BEC since its size char¬ 
acterized by <r(t) rolls down towards 0 as it tries to lower 
its potential energy. From the values provided in eq.(21), 
the onset of instability starts at N cr = 2.4 x 10 19 . While, 
if i?*(l) is smaller ss 1/im, the critical occupation number 
for 133 Cs rises to « 5 x 10 21 . However, when N < N cr , 
there is one positive root of eq.(22), and the correspond¬ 
ing excitation frequency is 2wo , albeit independent of 
n = 1 gravity. 


IV. CONCLUSION 


Within the ambit of variational method, we have found 
that occupation numbers in excess of N cr cause collapse 
of BEC for attractive gravity falling off as r~ 3 . This can 
be subjected to experimental verification only when one 
attains BECs with macroscopically large occupational 
numbers ss 10 19 - 10 22 . For higher values of N , even 
Newtonian gravity appears to have significant effect on 
the BEC dynamics that needs to be studied more care¬ 
fully. The consequences of n > 2 LED theories on BEC, 
though not covered in this paper, need to be studied. In 
particular, it would be interesting to see whether their ef¬ 
fects could be disentangled from those arising from other 
atomic interactions like van der Waals force. 
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